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We present the linear response matrix for a sliding domain wall in a rotatable magnetic nanowire, 
which is driven out of equilibrium by temperature and voltage bias, mechanical torque, and magnetic 
field. An expression for heat-current induced domain wall motion is derived. Application of On- 
sager's reciprocity relation leads to a unified description of the Barnett and Einstein-de Haas effects 
as well as spin-dependent thermoelectric properties. We envisage various heat pumps and engines, 
such as coolers driven by magnetic fields or mechanical rotation as well as nanoscale motors that 
convert temperature gradients into useful work. All parameters (with the exception of mechanical 
friction) can be computed microscopically by the scattering theory of transport. 

PACS numbers: 75.78.Fg,85.85.-|-j,62.25.-g,72.15. Jf 

I. INTRODUCTION 

Onsager's reciprocal relationJ^ reveal that seemingly unrelated phenomena can be expressions of identical micro- 
scopic correlations between thermodynamic variables of a given system.^^lThe archetypal example is the Onsager-Kelvin 
identity of thermopower and Peltier cooling. 

Research into the interaction between electric currents and the ferromagnetic order parameter o f th e last years 
has paid off handsomely. On one hand, the predicted charge current-induced spin transfer torqu^^ in metallic 
ferromagnetic structures such as spin valves and domain walla^"^ in ferromagnetic wires has been understood in some 
detail and applied to random access magnetic memories,'^ logics,'^ and shift registers.E^l On the other hand, it has 
been established that a moving magnetization pumps a spin currenlPi] that can be converted into charge currents and 
voltages in ferromagnetj normal metal bilayers,'^ ferromagnetic textures,'^ multilayerSj^i^ or by th e spin Hall effect.'^ 
Spin-pumping induced voltages have been observed in metallic magnetic heterostructures j^^^^^^ tunnel junctionJi^ 
and magnetic wires with a moving domain walL^i^ A modern illustration of the power of Onsager's relations is the 
demonstration that spin-transfer torques and charge currents induced by magnetization dynamics are two sides of the 
same medal .'^SHUl 

Domain walls also react to thermal gradients, as first observed and discussed by Jen and Berger.l^^HH Domain wall 
displacement by laser heating is a possible technology for high-density magnetic recording.'^ Hatami et alW^ and 
Saslo"w^ proposed a thermoelectric spin-transfer torque as mechanism for magnetization switching in spin valves and 
domain wall motion in magnetic wires. More recently, Kovalev et addressed this issue for general one-dimensional 
spin textures in ferromagnetic wires. Yuan et al. found large non-adiabatic corrections to the thermal torques on 
narrow domain walls .^^Sl 

The scattering theory of electron transport can be employed to describe dissipative processes in magnetic systems 
such as the Gilbert damping of magnetizations dynamics, ^'^ leading to a microscopic for malism for the Onsager 
coefficients that govern the interaction between charge currents and magnetization dynamics .1^312] 

Nearly a century ago it was discovered that in macroscopic bodies the magnetization of a ferromagnet couples to 
the mechanical degree of freedom: Barnett demonstrated that a mechanical rotation of a demagnetized ferromagnet 
creates a net magnetization along the rotation axis, ^'^ whereas Einstein and de Haas showed that reversing the magnetic 
moment of a ferromagnetic cylinder induces a mechanical torque.^"* Both effects are governed by the same gyromagnetic 
tensor.l^Sl The microscopic theory of mechanical and magnetic angular momentum coupling in nanostructures has 
recently been picked up again.l^SHall Moreover, Wallis et al^^ succeeded in measuring an Einstein-de Haas effect by 
agitating a magnetic cantilever. Zolfagharkhani et aim detected the mechanical torque induced by the decay of a 
current-induced magnetization, which can be interpreted as a variation of the Einstein-de Haas effect.SSlThe conditions 
to observe the Barnett effect in nanostructures have been estimated by Bretzel et al^^ 

In this paper we show that effects of magnetic, electric, thermal, and mechanical forces can be unified in a linear 
response matrix relating the conjugated thermodynamic variables for charge, energy, magnetization, and mechanical 
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rotation. In order to keep mathematics simple, we focus on a thin wire of an easy-plane metallic ferromagnet as 
studied by Hals et al^^ for current-induced magnetization dynamics. The wire is connects two heat and particle 
baths and is allowed to rotate (Figure 1). We may profit from Onsager's relations according to which we have to fill 
in only one half of the nondiagonal elements of the response matrix. This implies, for example, that in the linear 
regime Barnett and Einstein-de Haas effects are equivalent. We identify the heat-current driven domain wall motion 
and conclude that domain wall motion is associated with the pumping of heat (or a "thermal motive force"). The 
mechanical torque generated by a temperature difference opens the vista of magnetic nanoscale heat engines. 

In this paper we first recapitulate the basic thermodynamics following de Groot^ for a conventional thermoelectric 
element in Section Ila. In Section lib we show that the Onsager principles can be applied to the coupling between 
magnetic and mechanical dynamics for a model system of a magnetic wire containing a domain wall. In Section III the 
Onsager matrix is derived for a coupled thermoelectric and magnetomechanical system. In Section IV we specify how 
the Onsager matrix elements can be computed microscopically. Section V is devoted to a discussion of the magnitude 
of the couplings for a model system consisting of a nanowire of a ferromagnetic metal wire encapsulated in multO-wall 
carbon nanotubes. Section VI summarizes the conclusions. 



II. NON-EQUILIBRIUM THERMODYNAMICS 



According to the second law of thermodynamics the entropy S is maximal for the equilibrium state such that 

^ n n 

AS = - -^'^ gikaittk < (1) 
1=1 k=i 

for small deviations of the n state variables Ui = Ai — Ai from their equilibrium values Ai . The matrix of coefficients g 
is positive definite and symmetric. If we parameterize a small deviation of the system from thermodynamic equilibrium 
by the forces (affinities) Xi defined by (where T is the equilibrium temperature): 

X, = T^^-Ty^g,kak, (2) 

then, in linear response, the variables Ai will relax to their equilibrium values Ai according to 

n 

Ji = di — ^ LikXk, (3) 

k=l 

defining the response matrix L. Its elements can be introduced phcnomenologically or computed from microscopic 
principles by the Kubo- Greenwood formalism or scattering theory. The system responses Ji are called fluxes, currents, 
rates, velocities etc. Eq. ([3| remains valid in the presence of external forces slowly varying in time that may render 
the Ai time dependent. The entropy generation rate reads 

S ^ -^cii^gikak ^ Tj^^JiX^. (4) 

i k i 

Onsager discovered that, due to microscopic time-reversal symmetry, the linear response coefficients obey the reci- 
procity relations 

Lik (Hext, m) = eiekLkt (-Hext, -m) . (5) 

where = 1 if the state variable ai is even under time reversal and = —1 otherwise. Time-reversal (anti) symmetry 
in the presence of external magnetic fields Hcxt and equilibrium magnetic ordering indicated by a vector field with 
unit length m(r) (parametrizing the position-dependent direction of the magnetization) has been made explicit. The 
inverse of the response matrix L : 



X^^Y.^^k■h, (6) 



k=l 

has the same Onsager symmetry 

L'^k (Hcxt, m) = ^i^kL^l (-Hext, -m) • (7) 
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A. Thermoelectric element 



Consider as an example an ordinary thermoelectric element (such as a wire) connecting two reservoirs which are in 
respective thermal equilibria but at different temperatures T1/T2 and voltages Vi jVi. Let us define AT = Ti — T\ <C T 
and Ay = Vi — V\. If the wire has no independent degrees of freedom, we can describe a general (slightly out-of- 
equilibrium) state of this closed system by (half of) the energy and charge differences between the two reservoirs, 
\J = {U2 — Ui)/2 and q = ((72 — 9i)/2, respectively. Disregarding the wire's heat capacity and electrostatic capacitance 
relative to those of the large reservoirs, U and q correspond to the energy and charge that have been transferred from 
reservoir 1 (left) to reservoir 2 (right) with respect to some reference state. Jc — q and Jq — U are, respectively, charge 
and energy currents associated with U and q that are driven by AT and AV^. We next employ the thermodynamic 
identity 

(8) 



which holds for each reservoir separately. To leading order in the perturbations, the total entropy change S = Si+ S2 
introduced by moving a small amount of energy and charge between the reservoirs is thus 
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By comparison with Eq. Q, we identify the conjugate fluxes and forces: 

AT 

Jq = U , Xq -- 



T 



Jc^ q, Xc 



-AV. 



such that Eq. ^ becomes 
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The Onsager matrix can be rewritten in terms of the electric conductance 



G = 



heat conductance 



and thermopower or Seebeck coefficient 



such that 
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(9) 



(10) 
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Traditionally, the role of currents and voltages in the thermoelectric response are exchanged. In terms of the resistance 
R=l/G: 

' -AV' 
Jq 

Hereby, we recovered the Onsager-Kelvin relation between thermopower and Peltier coefficient: 

Jq 



R -TS 
TS Tk 



Jc 
AT 
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(15) 



(16) 



n = ^ =TS. (17) 

Jc AT=0 

The Sommerfeld approximation leads to the Wiedemann- Franz law k = GCT and Mott's formula S — 
—eCTdE^'nG\Ep, where Oe^t^GIep is the logarithmic energy derivative of the conductance at the Fermi energy 
Ep, C = (fcB/e)^7r^/3 is the Lorenz number and — e the electron charge. The dimensionless expression S'^/C vanishes 
quadratically at low temperatures, is small for most metals at room temperature,^^ and may usually be disregarded. 



Eq. ( 15 ) then becomes 
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1 TS 
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(18) 
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FIG. 1: Magnetic nanowire of length I in electrical and thermal contact with reservoirs. A domain wall is centered at position 
Tw The wire is mounted such that it can rotate around the a;-axis. A magnetic field and mechanical torque can be applied 
along X. 



B. Magnetomechanical element 

We consider now a quasi-one-dimensional magnetic nanowire with easy-plane anisotropy that contains a transverse 
domain wall, which is the standard model system for the study of magnetic domain wall motion. We chose here the 
tail-to-tail (rather than head-to- head) topology shown in Figure 1. The wire is mounted in a low- friction bearing 
such that it can freely rotate around its (a;) axis and a mechanical torque t™°'^'' can be applied. The system can 
also be driven by an applied magnetic field -ffcxt, and, via electric and thermal contacts, by a voltage (Ay) and/or 
temperature (AT) bias. 

Let us suppose initially that the magnetomechanical properties are decoupled from the electric and heat currents. 
The equation of motion of the magnetization MgXa {x,t) , where Af^ is the constant saturation magnetization, is 
governed by the Landau-Lifshitz-Gilbert (LLG) equation, appended by Barnett's gauge field that represents the 
aligning torque felt by angular momenta in rotating systems. In the frame of reference that rotates with the wireP^ 

rh = —7m X Hcff -I- am x rh + m x xip , (19) 

where 7 is the minus the gyromagnetic ratio (7 > for electrons) and ip the angular velocity of the wire around its 
axis. The effective field Heg is the functional derivative of the free energy J- with respect to the magnetization at 
rest, which has contributions from the applied, anisotropy, and exchange magnetic fields: 

SJ' [m] 

Heff = -YT^^-\ = (^cxt + Km^)ii - K^m.z + A^V^m , (20) 
Msom (r) 

where m = {mxjmyjmz) = (cos0,sin0cos0, sin^sint/)), the anisotropy constants K > Q, K± > 0, and the exchange 
stiffness A^^ have been introduced. In the absence of pinning the Walker ansatd^ 

In tan = - ^^^Zl^ and {x, t) = (/)(i) , (21) 

2 A^) 

provides a solution in terms of a domain wall with time-dependent position and (squared) width — Ac^l{K + 
K± sin^ (j)) . The polar angle </> is the tilt of the magnetization against the easy-plane anisotropy K± , which vanishes 
at equilibrium, (f) is a constant for sufficiently small, steady-state driving forces, and therefore not a dynamic variable 



in the regimes considered henceforth. Substituting Eq. (21) into the Landau-Lifshitz-Gilbert Eq. (19) 



, . IT \ ■ r,l 2((y9 - 7_ffoxt) 

rw = — (<^-7i/cxt), if_Lsm20 = . (22) 

a Q!7 

This solution is valid up to a critical (Walker) threshold field at which | sin 2^14/ 1 = 1. To linear order in the driving 
field, we can approximate the domain wall width A^ by its equilibrium value, A^; = \/ Aa^/K. 
The mechanical rotation of the wire is governed by the damped oscillator equation: 

/(^ + /3mcch^ ^ ^mcch^ (23) 
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where /3™<^<=*i is the mechanical damping parameter and r™'^'^'' the total mechanical torque acting along the a;-axis. The 
total angular momentum £axis of the mechanical and magnetic subsystems in a freely rotating wire of cross section A 

AM 

£axis = ^ (/ - 2r^ ) + /(p , (24) 

7 

is dissipated into the environment at a rate £axis = — This leads to an expression for the Einstein-de Haas 
torque induced by a moving domain wall: 



'EdH — 
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(25) 



We assume in the following that the system is overdamped, i.e., we limit our attention to frequencies smaller than 
PI I, such that the acceleration ip and moment of inertia / drop out of the problem. The rotation velocity is then 
directly proportional to the total torque r™°'^'^ = T™j?^^ + Tgjfg^: 

/3--V = r-r'^-— (26) 
7 

The mechanical energy E{ip) governs the external torque, t^^J?'^^ = —d^E{ip). 

The above results will now be shown to be consistent with Onsager's reciprocity principle and the second law of 
thermodynamics. Disregarding thermal effects, it is natural to switch to the free energy T instead of the entropy S: 

T{r^, ip)^T^+T^^ (2r„ - l)AMsH,^t + E{ip) , (27) 

where I is the total length of the wire and the domain wall position r^, is measured with respect to the left end of the 
wire. We omit the internal energy of the domain wall, which below the Walker threshold may be treated as a rigid 
particle-like mass-less object specified by its position. The conjugate forces associated with and ip are immediately 
found as 

X^ = ~ = -2AM,ffext ; = - = tZ""^ . (28) 



After simple algebra using Eqs. (22 1 and (26), the energy dissipation is found to be positive definite: 

TS = -f= -2AAhH,^,r^ + r^^V = + > , (29) 



Rewriting the equations of motion ( 22 1 and ( 26 ) , the cross terms are seen to obey Onsager's symmetry: 




1 + "^^1 wv) (30) 

The antisymmetry of the off-diagonal terms stems from Onsager's reciprocity, which relates here the response of the 
tail-to-tail domain wall to that of its time-reversed partner, which is a head-to-head domain wall. Note that the 



inverse of Eq. ( 30 1 is simpler 




mech 2AM 

2AM 2aIi a ) ( ^ ) 7 (31) 
7 7 A 



but Eq. (30) should be closer to experimental set-ups in practice. We may rewrite it as 




(32) 



where 

^mcch _ ^mcch _^ iX^AMg 

2KAM, 

a = a-\ T— . 

/Jmech 



(33) 
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The magnetomechanical coupling creates an apparently increased damping of the magnetization dynamics and/or the 
mechanical motion that is proportional to the number of spins in the domain wall. ^3™°'^'^ is the mechanical friction: 
when it becomes large the mechanical motion is quenched and the excess Gilbert damping is suppressed a — > a. In 
turn, the direct coupling of the mechanical torque to the rotation, L^^p = 1//?™'''^^, vanishes with vanishing Gilbert 
damping a, i.e., 7"™°"^'' is fully transferred into the magnetic system. For vanishing mechanical damping -|-j^g 
domain wall remains immobile under a magnetic field, but the wire rotates with an a-independent angular velocity, 
which exactly compensates the external field in the rotating frame (i.e., ip = ^Hcy±)- These results are valid only in 
the steady-state, overdamped mechanical regime considered here. 



III. MAGNETOMECHANOTHERMOELECTRIC SYSTEMS 



We now define the conjugate thermodynamical variables that allow us to take advantage of Onsager's relations 
as energy transfer U — {U2 — Ui)/2, charge transfer q — {q2 — 9i)/2, domain wall position r^,, and lattice-rotation 
angle ip. The corresponding fluxes are given by their time derivatives Jq = U , Jc = q, Jw = rw, Jip = ^- The 
thermodynamic forces ([2| depend in principle on the values of all thermodynamic variables. It is possible to work 
out a general scheme that includes all possible cross correlations, but it would not be very transparent. Instead, 
we follow a more pragmatic approach that is based on the low-temperature free energy for the magnetomechanical 
degrees of freedom, which are coupled to thermoelectric transport between the reservoirs by the spin torques. The 
linear response matrix then reads J = iX, where 

J = ( Jc, Jq, r-u. )^ (34) 
X = ( -AV, r-^^ -2AM,i/ext f (35) 

and 



L=\ ,'3"' . (36) 



Lcc 


LcQ 
























LwQ 







When thermoelectric and magnetomechanical systems are uncoupled, the matrix elements derived in Sec. II may 
be filled in unmodified. According to the Onsager symmetry, ^^^^(m) = Lu)a;(— m) = itoa:(m) and Lxip{va) — 
L^xi^'^) = ^Lipx{'va), for x = (c, Q), assuming that our system obeys a structural mirror symmetry with respect to 
a plane normal to the wire in Fig. 1. It is useful to introduce also the inverse matrix X = i^^J, recalling that L^^ 
and L have the same Onsager symmetry. 

We can draw a number of conclusions from the Onsager relations already, i^.c and Lew represent the Onsager equiv- 
alent pair of current-induced transfer torque and charge pumping by the magnetization dynamics, respectively.'^^"^ 
We know that a temperature gradient can induce a spin-transfer torque,!^ which is here represented by L^q- According 
to Onsager symmetry an opposite and equivalent eff'ect exists, i.e., a heat current induced by magnetization dynamics, 
which might be applied for cooling or heating purposes. As explained above, the mechanical motion induced by the 
magnetic field as quantified by L^p^ (Einstcin-de Haas effect) is identical with the Barnett response function —L^ip, 
which describes the magnetization dynamics induced by rotation (Barnett effect). S ince Ltpc — — L^^p the magnetic 
wire can be employed as an electromotoi^^ and electric generator. A temperature gradient induces a rotation of the 
wire via L^q, which leads to the prediction of a heat engine that can carry out mechanical work under a temperature 
difference. The opposite effect, in which mechanical motion of the wire is transformed into a temperature gradient is 
governed by Lq^ = -L^q. 

The remaining task is to work out the elements of the 4x4 response matrix. In the adiabatic regime, the magnetic 
texture varies slowly with respect to the magnetic coherence length 1/Ac = IV-^f ~ I/'^fL where X^p^'^ are the spin- 
dependent Fermi wavelengths and t / 4- denote the majority/minority spin carriers, respectively. The spin torque on, 
or angular momentum transfer to, the magnetization induced by an a pplie d voltage (superscript (0) indicates a static 
magnetization texture and the absence of thermoelectric effects) read^^^^^^^^ 

[T^^f^ = -^^PGAV{1 - /Sr^mx) ^m. (37) 



in terms of the spin polarization P = (G^ — G^) /G of the electric conductance G = G-^ -f- G^ of the single-domain 
ferromagnet. The prefactor (torque in the plane of the domain wall magnetization) can be easily identified as the 
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angular momentum rate of change of the spin-polarized carriers that corresponds to an adiabatic spin reversal. 
The out-of-plane torque components is caused by the mistracking of the spin in the magnetization texture that is 
parameterized by (3^'^^. The torque by the thermoelectric spin current induced by a temperature bias (superscript 
(0) again denoting a static texture and the condition S'^/C ^ 1) reads analogously: 

where P' = Oe (PG) /dsG is the polarization of the energy derivative of the conductance!^ The parameter (Sq^^ 
parameterizing the out-of-plane torque component differs from since the non-equilibrium energy distribution 

defining the spin current through the texture has a node at the Fermi energy rather than a maximum. The Seebeck 
coefficient for the homogeneous ferromagnet S = {S^G^ + SiGi) /G. 

In the coupled system, torques are induced by the magnetization and mechanical motion as well, which can be fully 
included into the equations of motion by adapting the charge and heat currents, rather than voltage and temperature 
as system variables (forces). The transformation can be carried out very generally, but leads to different parameters 
for the out-of-plane torques.l^By inverting the thermoelectric matrix in the Sommerfeld approximation and assuming 
S'^ / C ^ 1, the j3 parameters remain unmodified for the current biased torques, however. The torque induced by a 
heat current Jq then reads 

h ^ p' Q / \ f) 

mag / T I /:)mag ^ \ ^ 



e 2AMs CT 

whereas the charge current torque becomes 



jQ(l-/?-^W)£m, (39) 



= (-^e - P,) (1 - /^-mx) (40) 

Note that the conventional thermoelectric charge current has been subtracted here from the total charge current. 



Adding the spin torques (39) and (40 1 to the right-hand side of the LLG Eq. (19), we can employ the Walker ansatz 



again to solve for the current driven domain-wall velocity: 

A negative charge, thus positive particle, current and Pp^'^^ > pushes the domain wall to the right. For an 
electron-like thermopower {S < 0) and P' (3q^^ > 0, a positive heat current has the same effect. 



In order to relate (p to the mechanical torque and identify the unknown response coefficients in Eq. (36), we need to 



generalize the conservation of angular momentum, Eq. (26 ), to account for the spin currents injected into and drained 
from the leads by 



with, for the present choice of magnetization texture. 



^ST^' = -J(7cP (jc - +1q^Jq)- (43) 

The effect is maximized (7c = 7Q = 1) when the angular momentum is drained completely from the wire into the 
reservoirs. In the opposite limit, the spin currents are dissipated completely in the wire (rather than in the reservoir), 
as in sufficiently long normal metal terminals to the ferromagnet (N|F|N) that are part of the mounted wire, such 
that T'"j°'^'^ = 0. The domain wall equation of motion, including the Barnett torques induced by rotation, reads 



= fi?ext - + r„|j.,j« . (44) 
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From Eqs. ( 41|42|43|44[ ) we can specify all coefficients of the inverse response matrix 



L- 



. h pqma^ h S f p/ gmag pomagA 2AM s 2 AM s a 



(45) 



This representation appears to be less convenient for comparison with practical experiment. In a purely electric 
circuit it is possible to freely change from a current-biased to a voltage bias set-up. This appears less convenient for 
the other sets of conjugate variables. It is therefore necessary to adopt the results to the experimental problem at 
hand. For a set-up in which the driving forces are the Xi considered here, it is appropriate to invert the above matrix 
in order to obtain experimentally more relevant response functions. We have seen in the previous section that for 
the purely magnetomechanical system the inversion is equivalent to a renormalization of the damping constants. The 
inversion of the 4x4 matrix leads to lengthy expressions that cannot be interpreted that easily. The simplest approach 
is second order perturbation theory to estimate the importance of the self-consistent couplings. The diagonal elements 
of the response matrix then read 

L - ' ^1 (46) 
^" - 77^ P + 2- (T-i\ (T.-i\ ^ ' ^^^> 



while the non-diagonal elements become 



(47) 



In Eq. ( 45 1 the 2x2 thermoelectric matrix and the mechanical diagonal elements scale with the system length and 



inversely with the wire cross section I /A, whereas all others are independent of the system size. The non-diagonal 
block matrices may therefore be treated by perturbation theory in the long and/or narrow wire limit. By defining the 
block-diagonal matrix 










^incch 






2AMs 



\ 











2AMs 2AMs a 
7 7 A 



(48) 



and treating 5L ^ = L ^ — Lq as a perturbation, we find to lowest order in 6L 



Using the Sommerfeld approximation and letting S'^ /C 
as: 



Lk.Lq- USL-^U. (49) 
0, we obtain the elements of the lower non-diagonal block 



1 



2AM, a 



2AM, A„ 



7 P 



mcch 



(50) 



n 7 GST 
e 2AM, a 



jmag 



2AMs \w 
5moch 



7 /3" 



(51) 



h G 



mcch 



P[lc- 



(52) 



h GST . / PT' 



(53) 



The perturbation expansion holds well for the example treated below. If it turns out inaccurate, the full matrix should 
be diagonalized, of course. 

In strongly spin-orbit coupled systems the transport polarizations P,P' are not well defined.'^ However, the dynam- 
ics is governed only by the combinations PP^^^^, P'jc, F" Pq^^, P'jq, which can still be determined. ^'^ An important 
handle to experimental access to the different parameters is a comparison of heat-current driven domain wall motion 
for both closed and open electric circuits. In the limit i/^mcch _ jj^^^ = g, the former case leads to 
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whereas 

^ ~2AM—e a ) ^^^^ 

An interesting simplified system consists of a completely pinned magnetic domain wall. In this regime, the magnetic 
degrees of freedom drop out of the problem.-^ In the adiabatic limit, the spin current then transfers all angular 
momentum directly to the lattice. Since the magnetization does not move, there is no magnetic dissipation. The 
response functions in that limit (indicated by a prime) are obtained in the limit a — >■ cxd and are significantly simplified 



Kc^~^p^P^c, (56) 



Tf ^ GST . /~—\ 

L^Q = --^P^Q- (57) 



The Onsager equivalent to the current-induced rotation is the rotation-induced charge and heat pumping by the 
otherwise fixed magnetization texture in the domain wall. It can be explained in terms of the magnetization texture 
that carries out a rotation rather than a translation, which in the rotating frame results in an effective (Barnett-like) 
field 1^/7 between the two reservoirs, which drives the charge and heat currents. Whether the pinned or the moving 
magnetization more effectively transfer angular momentum between currents and lattice depends strongly on the ratio 
of the dissipative out-of-plane torques and the Gilbert damping constant. The situation in real domain walls with 
weak pinning will be somewhere between the extremes of rigid translation and full pinning, but its full treatment is 
beyond the scope of the present work. 



The dynamics of insulating ferromagnets can be obtained by simply crossing out the first row and column of Eq. ( 45 ) 
related to the charge degree of freedom. In the remaining 3x3 matrix, the spin torques are exerted by the pure 
heat currents carried by spin waves, unlike in metallic systems, in which the spin torque is dominated by the electric 
current. The detailed response function for insulating ferromagnets will be discussed separately, however. 



IV. SCATTERING THEORY 



The magnetic damping and the charge-current magnetization coupling have been determined microscopically by 
scattering theory.[22MllHere we briefly review the relevant published results and add new ones related to heat transport. 

The Onsager response functions derived above contain a number of parameters, basically the spin-dependent con- 
ductances at the Fermi energy = Go- {Ep) , the Gilbert damping a and the dissipative out-of-plane spin transfer 
torque associated to charge current /S™'*^ ^nd heat current /3q^^. They can all be written in terms of the scattering 

matrix S of the wire at a given energy.'^Slsil Using the conventional notation in terms of transmission (t, i') and reflec- 
tion (f, f') matrices^*^ the scattering matrix in the space of the transport channels to and from the wire at an energy 
E and spin indices ct, a' reads: 



The spin-dependent conductance of the (single-domain) ferromagnet can be expressed by the Landauer-Biittiker 
formula: 

G. = ^Tr^ (59) 

where the trace indicates the sum over (orbital) transport channels at the Fermi energy. The Mott formula for the 
conventional thermopower S = —eCTds In G can be computed from the energy-dependent conductance G = G^ + G^. 

When S is the scattering matrix of the ferromagnet including one domain wall at r^, , the parametric pumping of a 
charge current}^ by the moving domain walPl' into the right leads reads 

Jc,w = -^r^lTaTis-^S'^fz, (60) 
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where 

in the same space as the scattering matrix, i.e. propagating states in the left lead and the right lead, respectively, 
and Trs is the sum over these states (including spin). The expression for the energy pumped out of the system with 
a parametric time dependence of the scattering matrijP^HH 

h dsds^ h . ds ds^ 

has been employed by Brataas et aZpSl to derive microscopic expressions for the Gilbert damping and by Hals et al^ 
for the charge current-induced domain wall motion. When evaluating the scattering matrix for zero bias and assuming 
that the domain wall is driven by a magnetic field Je ~ i'^w)'^ /L^w- 

For S'^ / C <C 1 and absence of rotation the response function reads to lowest order in the conductance 

^'"-~2AM:a^^Ve2AM:a) ' ^^^^ 



e 2AMs a 

For long wire lengths / we recover the result by Hals et al^^ for the Gilbert damping.- 



7 '^p^mag^ 



7^A,„ dS dS^ 

a — - — — — hm Ir^— — — — . (56) 
Qi:AMs /-i-oo or^ju or^ 



and the dissipative torque correction 



p^mag ^ ^ iniTr,|^5tf,. (67) 

in /-i-oo G or,,, 



The heat current pumped by the magnetization dynamics depends linearly on the frequency and amplitude of 
the pumping parameter and should not be confused with the energy current Je , Eq. ( 62 1 , which is to leading 
order quadratic in these quantities. This thermoelectric contribution to the pumping current can be obtained by a 
Sommerfeld expansion of the energy dependent parametric pumping current as derived by Moskalet and Biittiker.'^ 
The heat current driven by a moving domain wall then reads 

jQ = -eCT^-^J,^^{E). (68) 



where Jc,w is a function of energy. Observing that the domain wall velocity in Eq. ( 60 ) is a parameter that can be 
pulled in front of the energy-derivative we arrive at 

This leads to 

lim L^Q = CT^ ^ — (70) 

In the limit of long wires the leading term of the heat-domain wall coupling 
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and therefore: 



P'PT' = .T"" , (72) 



which has been recently evaluated by Hals et al. for GaMnAs!^ We can also derive a relation between the parameters 
governing the charge current and heat current-induced domain wall motion 

where 

Sp = -e/:T^ In (GP/?,"'^s) . (74) 



NUMERICAL ESTIMATES 



To mount magnetic wires such that they can rotate freely seems challenging, but should be possible!^ Elias et al. 
have grown single-crystalline FeCo wires inside multi-wall carbon nanotubes. The outer walls of multi-wall carbon 
nanotubes form almost ideal bearings for the rotation of the inner tubes!^^'^ A possible recipe for creating a system 
that can be described by the present model is therefore a suspended bridge of a multi-wall coated FeCo nanowire. In 
order to insure that all currents flow through the ferromagnet, it might be useful to burn off the carbon in the free 
standing part. Such a system could sustain GHz rotation frequencies when driven by the spin-flip transf er t orque 
that dissipates an injected spin current The scaling with different material constant is obvious in Eq. (45 1. We 
chose parameters that are close to permalloy, viz. p = 10^^ 57 cm, A^, = 100 nm; S — — 40/iVK^^ . Servantie and 
GasparcPS] report the dynamic friction /3™°'^'^/Z = 0.044 unm/ps for a (4,4) nanotube rotating in a (9,9) nanotube 
bearing. We choose a wire area cross section A 100 nm^ and a wire length I = 1 /xm . We chose a damping of 
a — 0.01 and = P' [3q^^ = P^c^^^ = P'Iq^^ = 1- The response function becomes dimensionless by choosing 

appropriate units for the thermodynamic fluxes and forces. This leads to 



O.IT 



10-21 Nm 



O.IT 




-6.6 • 10- 

0.92 
-0.05 




W-'' J/s 

GHz 
10= m / s 



(75) 



or 



v 



W-'' J/s 

GHz 
10^ m / s 



AV 




(76) 



are 



We can make a number of observations. For the present example the self-consistency effects, e.g. L^^^ 7^ ^) 

well described by the perturbation approximation used earlier, since the off-diagonal block matrices coupling of the 
thermoelectric and magnetomechanical systems are rather small. These couplings can be increased by a large diameter 
or shorter length of the wire. In permalloy a temperature gradient of 0.2 Knm^^ induces a charge current density 
of lO^Acm"^, which should suffice to move the domain wall in state-of-the-art wires. A material with a smaller 
saturation magnetization and large dissipative torques such as GaMnAs will be more susceptible to heat and charge 
current-induced magnetization dynamics. The small friction of the nanotube-lubricated rotation causes the strong 
coupling between the mechanical degree of freedom and the magnetization dynamics. The best way to enhance the 
coupled dynamics is the use of materials with a low Gilbert damping, however. 



VI. SUMMARY, EXTENSIONS, AND CONCLUSIONS 



We derived the linear response matrix for a magnetic wire in contact with electric and thermal reservoirs that can 
rotate along its axis. 
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Jen and BergeJ^ observed domain wall motion in amorphous magnetic alloys under a temperature gradient as 
small as 0.1 K //Ltm from the hot to the cold side. They offer two alternative explanations, viz. an entropic driving 
force in a domain wall gasp2lor a domain wall drag by the eddy currents induced by the anomalous Nernst effect.l^In 
thin wires such as addressed here both mechanisms are unlikely to compete with the thermal spin-transfer torque.l^ 
The domain wall displacement due to temperature dependence of magnetic anisotropies as utilized by Miyakoshi et 
alW^ should not play a role in soft magnets such as permalloy, but temperature dependence of pinning potentials can 
affect the dynamics, in principle. 

Sinitsyn et al.^^ predicted a translational domain wall motion under rotation of the magnetization texture, finding an 
identical dependence of domain wall velocity with rotation frequency as we do.l^ However, Sinitsyn et alP^ consider 
damping in the laboratory frame of reference and not in the rotating frame. Their predictions therefore hold for 
domain walls rotated relative to the lattice with direct magnetic dissipation into the environment, whereas we focus 
here on combined rotations of lattice and magnetization with mechanical friction. 

We conclude that a moving domain wall pumps heat, which we might call domain wall Peltier effect. A sizable 
cooling power may be associated with magnetic-field induced domain wall motion. The domain wall drag by the 
thermal spin transfer as well as the domain wall cooling can be computed microscopically by the methods used 
by Hals et al.'^^ and Starikov et al^ Kovalev et al.^^ independently obtained results for the interaction between 
heat currents and magnetization. Proceeding from arbitrary one-dimensional textures they illustrate their results by 
a spin-spiral model rather than a single domain wall, however. Spin spirals should be sensitive to pinning effects 
due to (near) commensurability with the underlying lattice and at the wire terminals that are prone to suppress 
magnetization motion. For the heat and charge current driven magnetomechanical motion this can actually be an 
advantage. Kovalev et aJ^^ address thermal coolers by moving magnetization textures, but found low efhciencies for 
conventional magnetic materials. Hals et al. also point out that in GaMnAs the heating due to dissipation takes over 
any cooling effects already at moderate domain wall velocities. ^'^ 

The set-up in Figure 1 generates a charge current-induced mechanical torque by domain wall motion, which is 
quite different from the mechanical torque that is generated by a decaying spin accumulation,^!] or the spin-torque 
electromotor.'^ 

The spin-torque motors based on moving domain walls have a drawback: they can operate only with a single stroke, 
limited by the wire length over which the domain wall can propagate. A similar problem has been encountered for 
the DC electromotor, which has been solved by Faraday in the form of a commutator that periodically inverts the 
sign of the mechanical torque. However, a pinned texture (domain wall or spin spiral) as a rotor material solves this 
issue. Such a material would not profit from the enhanced out-of-plane dissipative torques predicted by Hals et al¥^ 

Many protein-based molecular motors in the cell may be Brownian motorJ^such as Feynman's ratchet and pawl,!^ 
in which stochastic motion in the presence of a temperature or chemical potential difference produces useful work. 
The present contraption also produces work out of a temperature difference on the nanoscale, thus can be interpreted 
as a realization of Feynman's ratchet, in which directionality provided by the "pawl" is replaced by the chirality of 
the ferromagnet. 

The present scheme can be extended into dilferent directions. An extension from one to two-dimensional textures 
is necessary to treat vortex domain walls in wider wires. The formalism is easily extended to describe the coupled 
motion of charges, lattice, energy and spins as a function of harmonic driving forces in the linear response regime. 
This would allow handling torsional vibrations that can be used to observe the basic phenomena more easily than a 
rotation.!^ When normal metal contacts are attached, spin currents and accumulations become explicit thermodynamic 
variables.!^ The spin-Seebeck effect!^ and its Onsager equivalent, the spin-Peltier effect, can then be handled. The 
Onsager relations in many-terminal structures such as those used in studies of the spin and anomalous Hall effects,!^ 
will be extended to the thermal counterparts, such as the spin and anomalous Nernst, Ettingshausen, and Righi-Le 
Due effects. 

In conclusion, we investigated the coupling of charge and heat currents with magnetization and lattice for a realistic 
model system. All parameters can be determined by independent experiments and are accessible to microscopic 
calculations. On the basis of the response matrix we predict various magnetic nanoscale heat engines and estimate 
the parameters that govern their efficiency. 

Acknowledgments 

We acknowledge helpful discussions with Alexey Kovalev and Kjetil Hals, who pointed out mistakes in the first 
version of the manuscript. This work was supported in part by the Dutch FOM Foundation, the Research Council 
of Norway, Grants Nos. 158518/143 and 158547/431, and EC Contract IST-033749 "DynaMax", the Alfred P. Sloan 
Foundation, and the NSF under Grant No. DMR-0840965. YT and GB would hke to thank the hospitality they 



13 



enjoyed at the TU Delft and NTNU Trondheim, respectively. 



1 L. Onsager, Phys. Rev. 37, 405 (1931) 

^ S.R. de Groot, Thermodynamics of irreversible processes (Interscience Publishers, 1952). 
^ J.C. Slonczewski, J. Magn. Magn. Mat. 159, LI (1996). 
* L. Berger, Phys. Rev. B 54, 9353 (1996). 

^ G. Tatara, H. Kohno, J. Shibata, Phys. Rep. 468, 213(2008). 

^ M. Klaui, J. Phys.: Condens. Matter 20 313001 (2008). 

'' D.C. Ralph and M. Stiles, J. Magn. Magn. Mater. 320, 1190 (2008). 

® R. Takemura, T. Kawahara, K. Miura, H. Yamamoto, J. Hayakawa, N. Matsuzaki, K. Ono, M. Yamanouchi, K. Ito, H. 
Takahashi, S. Ikeda, H. Hasegawa, H. Matsuoka, and H. Ohno, 2009 Symposium on VLSI Circuits, 16-18 June 2009, p. 
84-85. 

^ S. Matsunaga, J. Hayakawa, S. Ikeda, K. Miura, H. Hasegawa, T. Endoh, H. Ohno, T. Hanyu, Appl. Phys. Express 1, 091301 
(2008). 

^° S.S.P. Parkin, M. Hayashi, and L. Thomas, Science 320, 190 (2008) 

" Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, Phys. Rev. Lett. 88, 117601 (2002). 

^2 X. Wang, G.E.W. Bauer, B.J. van Wees, A. Brataas, and Y. Tserkovnyak, Phys. Rev. Lett. 97, 216602 (2006). 
" S.E. Barnes and S. Maekawa, Phys. Rev. Lett. 98, 246601 (2007). 

J. Xiao, G.E.W. Bauer, and A. Brataas, Phys. Rev. B 77, 224419 (2008). 

E. Saitoh, M. Ueda, H. Miyajima, and G. Tatara, Appl. Phys. Lett. 88, 182509 (2006). 

M. V. Costache, M. Sladkov, S. M. Watts, C.H. van der Wal, and B.J. van Wees, Phys. Rev. Lett. 97, 216603 (2006). 
^'^ T. Moriyama, R. Cao, X. Fan, G. Xuan, B. K. Nikolic, Y. Tserkovnyak, J. Kolodzey, and J. Q. Xiao, Phys. Rev. Lett. 100, 
067602 (2008). 

^® S.A. Yang, G.S.D. Beach, C. Knutson, D. Xiao, Q. Niu, M. Tsoi, and J.L. Erskine, Phys. Rev. Lett. 102, 067201 (2009). 

W.M. Saslow, Phys. Rev. B 76, 184434 (2007). 
^° R.A. Duine, Phys. Rev. B 77, 014409 (2008). 

21 Y. Tserkovnyak and M. Mecklenburg, Phys. Rev. B 77, 134407 (2008). 

22 L. Berger, J. Appl. Phys. 58, 450 (1985). 

S.U. Jen and L. Berger, J. Appl. Phys. 59, 1278 (1986). 
2* S.U. Jen and L. Berger, J. Appl. Phys. 59, 1285 (1986). 

2=5 T. Miyakoshi, Y. Miyaoka, T. Shiratori, Jap. J. Appl. Phys. Part 1, 43, 4906 (2004). 

28 M. Hatami, G.E.W. Bauer, Q. Zhang, and P.J. Kelly, Phys. Rev. Lett. 99, 066603 (2007). 

2'' A.A. Kovalev and Y. Tserkovnyak, Phys. Rev. B 80, 100408(R) (2009). 

2* Z. Yuan, S. Wang, and K. Xia, Sol. Stat. Commun., accepted for publication. 

2^ A. Brataas, Y. Tserkovnyak, and G.E.W. Bauer, Phys. Rev. Lett. 101, 037207 (2008). 

^° A. Starikov, P.J. Kelly, et al. (unpubhshed). 

^1 K.M.D. Hals, A.K. Nguyen, and A. Brataas, Phys. Rev. Lett. 102, 256601 (2009). 
^2 K. M. D. Hals, A. Brataas, Y. Tserkovnyak, arXiv:0905.4170 

33 S. J. Barnett, Phys. Rev. 6, 239 (1915); S. J. Barnett, Rev. Mod. Phys. 7, 129 (1935). 
3^ A. Einstein and W. J. de Haas, Deutsche Physikalische Gesellschaft, Verhandlungen 17, 152 (1915). 
3^ L. D. Landau, E. M. Lifshitz, and L.P. Pitaevskii, Electrodynamics of Continuous Media (Pergamon, Oxford, 1984). 
P. Fulde and S. Kettemann, Ann. Phys. 7, 214 (1998). 

A. A. Kovalev, G. E. W. Bauer, A. Brataas, Appl. Phys. Lett. 83, 1584 (2003). 
38 A. A. Kovalev, G. E. W. Bauer, A. Brataas, Phys. Rev. Lett. 94, 167201 (2005). 

A. A. Kovalev, G. E. W. Bauer, A. Brataas, Phys. Rev. B 75, 014430 (2007). 
■1° T. M. Wallis, J. Moreland, and P. Kabos, Appl. Phys. Lett. 89, 122502 (2006). 

G. Zolfagharkhani, A. Gaidarzhy, P. Degiovanni, S. Kettemann, P. Fulde, and P. Mohanty, Nat. Nanotech. 3, 720 (2008). 
■12 A. A. Kovalev, Nat. Nanotech. 3, 710 (2008). 

■13 S. Bretzel, G. E. W. Bauer, Y. Tserkovnyak, and A. Brataas, Appl. Phys. Lett. 95, 122504 (2009). 

N. L. Schryer and L. R. Walker, J. Appl. Phys. 45, 5406 (1974). 

S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 
■"^ A. Thiaville, Y. Nakatani, J. Miltat, and Y. Suzuki, Europhys. Lett. 69, 990 (2005). 
■'■^ A. Mougin, M. Cormier, J. P. Adam, P.J. Metaxas, and J. Ferre, Europhys. Lett. 78, 57007 (2007). 
■1* Y. Tserkovnyak, A. Brataas, and G.E.W. Bauer, J. Magn. Magn. Mat. 320, 1282 (2008). 

In the limit of strong spin orbit scattering in which spin currents and spin-dependent conductivities are not well defined, 

Eqs. (|39| and (40 1 still hold by intepreting the prefactors as leading terms in a formal gradient expansion of charge and heat 



Y. V. Nazarov and Y. M. Blanter, Quantum Transport, (Cambridge University Press, Cambridge, 2009). 
M. Biittiker, H. Thomas, and A. Pretre, Z. Phys. B 94, 133 (1994). 
^^2 J. E. Avron, A. Elgart, G. M. Graf, and L. Sadun, Phys. Rev. Lett. 87, 236601 (2001). 



14 



K. M. D. Hals, G.E.W. Bauer, and A. Brataas, Sol. Stat. Commun., in press. 

It should be noted that the accuracy of the Walker ansatz for driven but pinned domain walls assumed here, i.e. disre- 
garding the internal degree of freedom of the domain wall texture including its interaction with spin waves, requires further 
investigations. 

M. Moskalets and M. Biittiker, Phys. Rev. B 66, 035306 (2002). 

A. L. Elias, J. A. Rodriguez-Manzo, M. R. McCartney, et al., Nano Lett. 5, 467 (2005). 
A. M. Fennimore, T. D. Yuzvinsky, W. Q. Han, et al, Nature 424, 408 (2003). 
J. Servantie and P. Gaspard, Phys. Rev. Lett. 97, 186106 (2006). 

N.A. Sinitsyn, V.V. Dobrovitski, S. Urazhdin, and A. Saxena, Phys. Rev B 77, 212405 (2008). 
A. A. Kovalev and Y. Tsorkovnyak, Sol. Stat. Commun., in press. 

R. D. Astumian, Science 276, 917 (1997); P. Hanggi and F. Marchesoni, Rev. Mod. Phys. 81, 387 (2009). 

R. P. Feynman, The Feynman Lectures on Physics, Vol. 1, Ch. 46, ( Addison- Wesley.Massachusetts, 1963). 

K. Uchida, S. Takahashi, K. Harii, ,J. Icda, W. Koshibac, K. Ando, S. Mackawa, and E. Saitoh, Nature 455, 778 (2008). 

E. M. Hankiewicz, J. Li, T. Jungwirth, Q. Niu, S.-Q. Shen, and J. Sinova, Phys. Rev. B 72, 155305 (2005). 



